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Abstract. We first review the nonrelativistic lagrangian theory as a 
framework for the MOND equation. Obstructions to a relativistic ver- 
sion of it are discussed leading up to TeVeS, a relativistic tensor- vector- 
scalar field theory which displays both MOND and Newtonian limits. 
The whys for its particular structure are discussed and its achievements 
so far are summarized. 
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, The success of Milgrom's MOND paradigm (Milgrom 1983) in explaining the 
' shapes of disk galaxy rotation curves, and correlations between the parameters of 
^ r*| such systems (e.g. Tully-Fisher law) with no appeal to dark haloes is by now well 
O !' known — for a modern review see Sanders & McGaugh (2003; henceforth S&McG). 
q , This nonrelativistic scheme is summed up by a novel relation between acceleration 
!_i a of each element of an extragalactic system and the local Newtonian field — V<I>,/v, 

M(|a|/a )a= -V$ N , (1-1) 

,£h ! where ao ~ 10~ 10 m s~ 2 , and the \x function interpolates between \i(x) = x for 
■ x <C 1 and /i(oo) = 1. It is this last proviso which guarantees correspondence 
H \ of MOND with Newtonian gravity in the laboratory or the solar system, both of 
■ - - i which are supra ao situations. 

MOND lends itself to reformulation as a gravitational field theory alternative to 
Newton's. In this lagrangian based theory, named AQUAL (Bekenstein & Milgrom 
1984; henceforth B&M), the Poisson equation is replaced by 

V ■ bj(|V$|/a )V$] = AnGp, (1.2) 

with <!> the physical gravitational potential (so that a = — V$). In spherically 
symmetric situations Eq. (|1.2f> has Eq. (|l.lf> as a unique solution. But when 
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symmetry is lower, an extra term appears on the r.h.s. It is this corrected MOND 
equation which is to be used in place of Eq. deriving as it does from a 

lagrangian, it guarantees that momentum, energy and angular momentum will be 
conserved in MOND dynamics. The original MOND formula failed to assure this. 
AQUAL has a number of further advantages over plain MOND (B&M, S&McG). 

AQUAL is thus a fruitful tool for exploring non dark matter (DM) resolutions 
of the missing mass problems. But it has shortcomings. For example, like the 
naive MOND formula, it fails to do away with the need for DM in clusters of 
galaxies, although its requirements for the dark stuff are certainly lower than in 
Newtonian gravity. And since it is basically nonrelativistic, AQUAL is unable to 
address a number of key subjects with relativistic overtones: gravitational lensing, 
cosmological evolution and the growth of structure, to name a few. How to cast 
the MOND idea in a relativistic mold while retaining the benefits of AQUAL ? 

2 Relativistic MOND theories leading up to TeVeS 

A relativistic version of AQUAL (call it RAQUAL) was offered in its defining 
paper (B&M). The basic idea was to replace the physical potential $ by a scalar 
field 4> whose equation is a relativistic reflection of Eq. 11. 2|) . This equation is 
written on a metric g a /3 obeying Einstein's equations (with a contribution to the 
energy momentum tensor from (f). One then defines the physical metric as 

g aP = e 24, g a f3. (2.1) 

This last describes a spacetime which has been locally stretched in every direction 
with respect to that described by g a p. All matter and radiations are supposed to 
sense only g a fl, and the requisite equations must contain only g a p and never g a fj 
by itself. RAQUAL reproduces nonrelativistic AQUAL, and so recovers much of 
the successful MOND phenomenology. 

However, RAQUAL fails on two counts. It permits superluminal propagation 
of tfi waves (B&M). And it is unable to give an account of gravitational lensing in 
agreement with the basic observation that lensing by galaxy clusters is anomalously 
strong compared to what was to be expected in view of their galaxies and gas 
content. This last problem is rooted in the conformal relation H2.1JI between the 
Einstein and physical metrics. Conformally related metrics share the same tracks 
for photons (or electromagnetic wave packets). And because 0's contribution to 
the source of Einstein's equations is weak in clusters of galaxies, the metric relevant 
for lensing in RAQUAL is nearly that in GR with no DM. 

Breaking the conformal relation has been the central theme in attempts to fix 
this problem. For example, it is possible to add to Eq. 1|2.1|1 the symmetric tensor 
4>.a 4>,0 multiplied by some scalar (Bekenstein & Sanders 1994). It turns out that 
in order for gravitational waves — perturbations of g a f3 — to propagate causally, i.e. 
within the light cone of g a f3, the mentioned scalar factor must be negative, and it is 
then found to weaken gravitational lensing, rather than enhancing it as intended. 

A way out of the impasse (Sanders 1997) is to introduce a constant vector field 
which points in the time direction, and then to stretch spacetime along the vector 
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by a factor e 2 ^ while shrinking it by the same factor orthogonally to the vector. 
This does the trick: provided the scalar field (j> used comes from a RAQUAL-type 
equation, MOND phenomenology is recovered while the lensing is augmented to 
the right proportions. However, this "Stratified" theory is not covariant (one needs 
the vector to be in the time direction, and this cannot be true in all coordinate 
systems). Covariance of equations of a physical theory is so ingrained in modern 
physics that we have to look further. 

3 The TeVeS field theory 

TeVeS is a tensor-vector-scalar covariant field theory (Bekenstein 2004; henceforth 
B04) incorporating various motifs already mentioned. Just as in RAQUAL, the 
tensor is an Einstein metric g a @ out of which is built the standard Einstein- Hilbert 
action. One passes to the physical metric as in the stratified theory, but now using 
an everywhere timelike dynamical normalized vector field U a . In equations 

g a p = e- 2<t, (g aP +U a Up) - tP+UJAp, (3.1) 

together with the normalization requirement g a ^U a Up = — 1. 
The dynamics of U a is dictated by the action 

S v = J [g afl g^U M U [M - 2{X/K){g^U ll U v + 1)] (-g) 1/2 d 4 x, (3.2) 

where the square brackets denote antisymmetrization. The K here is a dimension- 
less coupling constant, one of the parameters of the theory. The A is a Lagrange 
multiplier field included to enforce the normalization of the vector field; A is to 
be determined as the equations are solved. The kinetic part of S v is that of an 
Abelian gauge field, but the unit norm constraint deprives lA a of the proverbial 
gauge freedom. In addition, U a is not coupled to a current of any sort, so it would 
be incorrect to view the vector as mediating a repulsive Coulomb-like force. 
The dynamical scalar field <fi required in Eq. Il.'l 111 is assigned the action 

S s = -\J [^ Q/ V,«</\/3 + ^G£- 2 a 4 F(kGa 2 )] (-g) 1/2 d 4 x, (3.3) 

which involves an auxiliary nondynamical scalar field er, a new dimensionless pa- 
rameter, fc, a scale parameter, £, and an unspecified function F. In the kinetic part 
of S s , <j),a4>,0 is contracted not with g a/3 , but with the tensor h a>3 = g a/3 — U a U^\ 
this is introduced to cope with the problem of supcrluminal propagation of scalar 
perturbations evidenced by RAQUAL. With h a ° the propagation is found to be 
subluminal with respect to the physical metric on condition that cf> > (B04). 
Contrary to appearances, S s is not a quadratic action. First a 2 must be expressed 
in terms of <f>^ a by varying S s with respect to a 2 , and inverting the resulting equa- 
tion. When the result is substituted in both terms of S s , the action becomes 
aquadratic in cj) a , resembling RAQUAL's action, with the scalar h a ^(f> :a (f> p re- 
placing g a ^4>^ a 4* (3- The form l|3.3fl is, however, more convenient for calculation. 
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These pieces of the total action must be supplemented with the matter action 
iS TO . In accordance with the Einstein equivalence principle (Will 1993) it must be 
formed by replacing g a p i— > g a p everywhere in the matter action in use in general 
relativity (GR), and modifying all covariant derivatives accordingly. In this way 
all types of matter respond to gravitational fields in the same way, and it could 
be said that matter and radiation delineate the physical metric. Now because g a p 
contains <\> as we U as Ma, the variations of the total action with respect to these 
quantities engender variations of S m which create sources for the equations of 
and U a . No direct coupling of these fields to the matter variables is thus required, 
which is pleasant as no principles are evident for determining such couplings. 

Is TeVeS unique ? Not at all. For example, one could choose to write S s 
and S v with g a p instead of g a p. The present approach seems preferable in that 
evidently g a p, 4> and lA a are all gravitational fields, so their actions should not 
involve the "matter's metric" g a /3- Further, one could add to S s a kinetic term 
for a such as is included in Phase Coupled Gravity (PCG), a theory devised in 
a partially successful attempt to forestall superluminal propagation (Bekenstein 
1986, 1990). In addition one could redefine h af3 as g a/3 - \u a U^ (C should be 
positive but not small in order to obviate superluminality; it could be a constant). 
An additional freedom is the choice of function F in Eq. 1)3. 3|) : each such choice 
defines a distinct theory. Some of the freedoms mentioned have been used to make 
a variant of TeVeS (Sanders 2005a, 2005b; henceforth S05). 

4 Results from TeVeS 

TeVeS has GR as a limit for K -> and I -> oo (B04). Thus if \K\ < 1 and the 
size of the system of interest falls well short of £, predictions of TeVeS and GR 
should not be that different. To verify that TeVeS can also give MOND in some 
other limit one must make a choice for the function F in action 1)3.3(1 . B04 makes 
the choice which is equivalent to the relation 

k£ 2 h^(/)^(j>^ = (3/4)fc 2 G*V(fcGcr 2 - 2) 2 (1 - kGa 2 )- 1 . (4.1) 

This is to be regarded as defining a toy theory. A better choice would have the 
r.h.s. change sign smoothly through kGa 2 = 1, a transition corresponding to the 
passage between a quasistatic system such as a cluster of galaxies {h^ u 4>^4> M > 0) 
and cosmology (h^ v ^) lfi cf) lU < 0) (see B04). 

Imagine linearizing the Einstein-like equations of TeVeS about flat spacctime 
(ga/3 = Vaf3 + • ' ■) hi a quasistatic situation. The final result for the physical metric 
after suitable choice of units is (B04) 

ds 2 = -{l + 2<$>)dt 2 + {l-2<$>)5 lJ dx l dx j ] $ = Z<P N + (f>. (4.2) 

Here $tv is the usual Newtonian potential of the baryonic matter and 5 is a con- 
stant depending on K and on the asymptotic (cosmological) value of the scalar 
field 4>; S is expected to be close to unity. With $ interpreted as the physical 
gravitational potential, metric l|4.2(l is precisely that used in GR to discuss both 
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dynamics and gravitational lensing. An immediate conclusion is that if GR suc- 
ceeds in describing both the dynamics of, say, a galaxy, and its lensing properties 
with a single assumed DM distribution, and if TeVeS succeeds in predicting the 
dynamics everywhere from the observed baryonic distribution, then TeVeS will 
predict the same lensing pattern as GR. Thus the problem that accosted RAQUAL 
and variants of it (Sec. is solved: lensing can be described using MOND ideas. 
This much was clear earlier from the non-covariant framework (Sanders 1997). Of 
course the problem with galaxy cluster dynamics (Sec. Q still looms. But it is 
important to realize that the two mentioned problems are distinct. 

Obviously, to determine $ requires one to calculate <p from the baryonic mass 
distribution using the TeVeS analog of Eq. 1)1. 2|) . In the linearized framework 
mentioned and assuming spherical symmetry, one can be immediately integrate 
the said equation once to get (we assume k -C 1 as well as 3 w 1) 

^V$ = V$at, /i« v a « —— . (4.3) 

1 + v/l + 4|V$|/a 4tt£ 

Thus we get the MOND equation as well as a formula for /i showing the 

appropriate asymptotic behavior at small and large |V$|/ao, and the value of do 
expressed in terms of the TeVeS parameters. MOND is recovered. Famaey and 
Binney (2005) have noted that the Galaxy's rotation curve is not fully consistent 
with this /i. Apparently a better F is needed. 




log r (A.U.) 



Fig. 1. |V$| (dashed), |V($ - $jv)| (in units of 1(T 8 cms" 2 ) with k = 0.03 (solid) or 
k = 0.1 (dotted) from TeVeS plotted vs the distance from the Sun. Observed constraints 
on the non-Newtonian part of the acceleration, |V(<& — $jv)|, are (left to right): from 
the precession of perihelion of Mercury, and of Icarus, from variation of Kepler's con- 
stant between Earth and Mars, between inner planets and Jupiter, Uranus or Neptune, 
respectively. The solid bar is the Pioneer anomaly range. 

The approximation inherent in Eq. i|4.3f) is good for all extragalactic applica- 
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tions, but begins to fail when |V$| > (7r/fc) 2 ao. In such situations, e.g. the solar 
system, it is best to start from the first term of the Laurent expansion of the r.h.s. 
of Eq. H4.1fl about kGa 2 = 1 . One then obtains 



The fact that /i asymptotes to (S+fc/47r) _1 as |V<&| — * oo means that for arbitrarily 
strong fields gravity is Newtonian, except that the effective gravitational constant 
is Gat = (S + k/4ir) G. It should be mentioned that the salient post-Newtonian 
parameters (Will 1993), important for inner solar system celestial dynamics, agree 
with those of GR; in particular (3 = 1 and 7=1 (Giannios 2004, B04). 

Notice the slow approach of \i to its asymptotic value as |V<I>| strengthens. 
One consequence is a slow increase of the Kepler constant a 3 /P 2 with distance 
from the Sun. Fig. 1 illustrates this effect and compares it with the constraints 
set by measurements of the orbits of diverse planets. Although a better form 
of F is clearly required, the significant point is that, in such theories, the total 
force in the outer Solar System will inevitably deviate from inverse square. This 
is of considerable interest in view of the Pioneer anomaly (Anderson et al. 1998, 
Turyshev 2005). 
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